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> ■ ABSTRACT 
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. Context. The observation of young stars with circumstellar disks suggests that the disks are dissipated, starting from the inner region, by the 
\l ' radiation of the central star and eventually by the formation of rocky planetesimals, over a time scale of several million years. It was also 
, shown that strong UV radiation emitted by nearby massive stars can heat a circumstellar disk up to some thousand degrees, inducing the 
photoevaporation of the gas. This process strongly reduces the dissipation time scale. 

Aims. We study whether there exists a correlation between the spatial distribution of stars with circumstellar disks and the position of massive 
stars with spectral class earlier than B5, in the open cluster NGC 661 1. 
' Methods. We created a multiband catalog of the cluster, down to V~23"', using optical data from a WFI observation at 2.2m of ESO in the B VI 
Or bands, the 2MASS public point source catalog and an archival X-ray observation made with CHANDRA/ACIS. We selected the stars with 
q , infrared excess (due to the emission of a circumstellar disk) using suitable color indices independent of extinction, and studied their spatial 
'_h ' distribution. 

rj^ Results. The spatial distribution of the stars with K band excess (due to the presence of a circumstellar disk) is and correlated with that of the 
■ massive stars: the disks are more frequent at large distances from these stars. We argue that this is in agreement with the hypothesis that the 
circumstellar disks are heated by the UV radiation from the massive stars and photoevaporated. 

•rH . 

£H Key words, stars: formation - planetary systems: protoplanetary disks - stars: pre-main sequence - open cluster and associations: individual 
. 5t , (NGC 6611) - infrared: stars - X-ray: stars 



1. Introduction 

Circumstellar disks are important features of star formation 
process for many reasons. For example, they regulate the an- 
gular momentum evolution of the stars and can evolve in plan- 
etary systems. Disks form during the gravitational contraction 
of the protostellar nebula and are dissip ated over time sca le of 
several millions years (see, for example, Hai sch et all2 001) . 

When it is not possible to directly observe the circum- 
stellar disks, as in the nearest star-forming regions, the disk 
dust emission in the infrared bands is used to detect and 
study them. However, observations show a decline in the emis- 
sion of the inner region of the disks (mostly NIR and mid- 
infrared ban ds) between ~1 and ~10 Myr ( see, for example, 
ISicilia-Aguilar et al.l20ollCalvet et alJ2005I) . This decline can 
be explained as the normal evolution of the disks, disrupted first 
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* Based on observations made with the European Observatory tele- 
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in the inner region by the radiation of the central star and by 
the formation of rocky planetesi mals, that are thought to occur 
fastest in the inner disks ( Podolak et al l ll993fc . Stars with tran- 
sition disks that are not characterized by NIR emission but that 
have strong mid- and far-infrared emission have been observed. 

The evolution of circumstellar disks, however, depends 
strongly on the environment where they formed. It was ob- 
served that strong UV radiation (mostly the far UV, with 
photons with energy between 6.0 eV and 13.6 eV) from 
the massive stars can heat the temperature of circumstellar 
disks up to thous ands of degrees, inducing the photoevapo- 
ration of the gas faollenbachl Il994 LTohnstone et all [l998; 

£ |A 



IStorzer & Hollenbachl Il999t lAdams et all 12004 . The disk 
mass loss rate is i ndependent of the distance from the cen- 
tral stars faollenbachl 1 1994 and since the surface de nsity de- 
crease s with increasing distance from the central star ( Pringlel 
1981), the disks are disrupted starting from the outer regions. 
So the evolution and the thermal structure of disks near massive 
stars is very different from those of isolated disks. For exam- 



2 



M. G. Guarcello et al.: Spatial distribution of the circumstellar disks in NGC 6611 



pie, Ijohnstone et alJ dl998h shown that the outer region of a 
circumstellar disk with a mass equal to 0.2 solar masses and a 
distance less than 10 17 cm from the star 1 Ori (spectral class 
06) is disrupted in 10 6 years up to a radius of about 1 AU from 
the central star. 

To investigate such effects, we study the stars with circum- 
stellar disks in the young open cluster NGC 6611. This clus- 
ter is characterized by a population of pre-main sequence stars 
younger than 3 million of years (see Sect. 14. 2t and by a large 
number of massive stars (56 star s with spectral class earlier 
than B5, Iffillenhrand etaiJll993l) . The massive stars are dis- 
tributed irregularly in the region, therefore the UV flux is non- 
uniform in the cluster. NGC 6611 is an ideal target to study 
the effects of UV radiation on the evolution of the circumstel- 
lar disks of low mass stars. The cluster is also within the Eagle 
Nebula that strongly absorbs the radiation of background stars. 
This will be very helpful to estimate the distance of the cluster 
and to avoid all the problems related to contaminating back- 
gro und stars. 

iHillenbrand etafl il993h (thereafter H93) estimated the 
age of the pre-main sequence stars of NGC 6611 as be- 
tween 0.25 and 1 million of years, the cluster distance as 
2.0 ± 0. 1 kpc and derived an anomalous extinction law, with 
R=3.75. iBelikov et alJ i 19991) (thereafter B99) derived a dis- 
tance of 2. 14 ± 0. 10 kpc, while smaller distance values are re- 
ported in the more recen t Catalog of Galactic Open Clusters of 
[ Kharchenko et alJ J2005h (1.7 kpc) and in the 2MASS study of 
Bonattoet alJ d2006h (1.8 + 0.5 kpc). 

We generated a multiband catalog using optical, near- 
infrared and X-ray data. Optical data will be used to determine 
the photospheric properties of the stars; NIR data will give indi- 
cations of the presence of circumstellar disks; X-ray data will 
be used as a membership criterion, particularly important for 
diskless cluster members. 

In Sect. 13 we will discuss the reduction of the optical data; 
in Sect. we will present NIR and X-ray data and the realiza- 
tion of the multiband catalog; in Sect. |4] we will use the op- 
tical Color-Magnitude diagrams to estimate distance, average 
extinction and age spread of the cluster; in Sect.[5]we will iden- 
tify all the stars with circumstellar disks; in Sect.[6]we defined 
the cluster members and studied the structure of NGC 66 1 1 
and in Sect. we will verify the effects of the UV radiation of 
the massive stars on the spatial distribution of the circumstellar 
disks. 

2. Optical data 

The optical data have been taken in the BVI bands with the 
Wield Field Camera (WFI), mounted on the 2.2 meter tele- 
scope of the European Southern Observatory (ESO) in La Silla 
(Chile). 




Fig. 1. WFI image in band I of NGC 661 1. The exposure time 
was 240 sec. The figure also shows the ACIS FOV, centered in 
the cluster, and the identifications of the CCDs of WFI. 



of each CCD is then 8'. 12 x 16' .25, and the total WFI FOV is 
33' x34'. 

The observations are part of the program ESO Imaging 
Survey (EIS), includ ing also the PRE-FLAMES program 
dMomanv et allkOOlh . Details of the 10 images of NGC 661 1 
are given in Tabled the seeing conditions were good during 
all the night, which was of photometric quality. Figure^shows 
a WFI image in I band, with the ACIS FOV (see Sect.Q also 
indicated. The cluster is in the center of the image. In addition 
19 images of 3 standard fields were used for the photometric 
calibration (see Sect. 12. 2\ . 

Instrumental calibration (Bias, Overscan and Flat-Field 
corrections) was performed with IRAF using the MSCRED pack- 
age. It was also necessary to remove the Fringe Pattern from 
the images in I band. We used the IRAF task RMVFRINGE and 
the WFI fringe pattern 1 that reproduces the interference figure. 

W e used t he DA0PH0T II/ALLSTAR/ALLFRAME proce- 
dures iStetsonill987l Il994l) to detect the point sources in the 
WFI images and to obtain their BVI magnitudes and their po- 
sitions from the PSF fitting. For all exposures, we also ob- 
tained an aper ture correc tion to the magnitudes with the routine 
DAOGROW ( Stetson, 1990) which uses the growth-curve method. 

2.2. Photometric calibration 



2.1. Optical observations and data reduction. 

The optical images used in this work were taken on 29 July 
2000 with WFI, a mosaic of 4 x 2 CCDs of 2048 x 4096 pixels, 
separated by gaps of 23". 8 in right ascension and 14". 3 in dec- 
lination, for a filling factor of 95.5%. The field of view (FOV) 



The instrumental magnitud es of stars in the standard fields 
SA107, SA92 and SA112 llLandoltl fl99l IStetsonl Eoool) were 
obtained using DAOPHOT II/ALLSTAR. 



1 available at: http://www.ls.eso.org/lasilla/sciops/2p2/E2p2M/WFI/ 
CalPlan/fringing 
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Image 




Exposure time(sec) 


Starting Time 


Filter 


Seeing(ar 


NGC 6611 


- 1 


30 


00:17:02.136 


Ic/lwp 


1.17 


NGC 6611 


-2 


240 


00:19:18.375 


Ic/lwp 


1.31 


NGC 6611 


-3 


240 


00:25:00.508 


Ic/lwp 


1.07 


NGC 6611 


-4 


30 


00:31:28.294 


V/S9 


1.43 


NGC 6611 


-5 


240 


00:33:43.779 


K/89 


1.31 


NGC 6611 


-6 


240 


00:39:21.946 


V/89 


1.43 


NGC 6611 


-7 


30 


00:46:09.205 


B/99 


1.43 


NGC 6611 


-8 


30 


01:00:40.413 


B/99 


1.07 


NGC 6611 


-9 


240 


01:02:21.695 


B/99 


1.07 


NGC 6611 


- 10 


240 


01:07:40.906 


B/99 


1.17 



The transformations from the instrumental magnitude sys- 
tem to the standard Johnson-Kron-Cousin one were computed 
for each CCD using the routine CCDSTD; the transformations 
were then applied to the standard stars by the routine CCDAVE 
and finally to all the stars by the routine TRIAL rtStetsonlEooll) . 

Usually, the typical form of the transformations for the pho- 
tometric calibration is: 

O-M + Z + KxQ + Ai xC (1) 

where O are the instrumental magnitudes of the stars; M are 
the standard magnitudes in the Johnson-Kron-Cousin system, 
Z is the magnitude zero-point, Q is the airmass of the observa- 
tions, C is a suitable color; Z, A\ and the extinction coefficient 
K are the coefficients to be determined. 

In our case it was necessary to include a position polyno- 
mial in the transformations, to allow for the dependence of 
the observed magnitudes on the position in the WFI mosaic, 
whic h is due to the non uniform illumination of the detec- 
tor jManfroid & Selmarl EoOll iKoch etall 120031 Icorsi et all 
2003$. Since the standard fields were not observed in a suf- 
ficiently wide interval of airmasses (from 1.148 to 1.22), it 
was also necessary to impose the extinction term's coefficients 
equal to the typical values for La Silla 2 : 0.14, 0.28 and 0.09 
respectively for the V, B and I band. We also included a tem- 
poral term in the transformations to take into account the time 
variability of the extinction. Therefore, the general form of the 
transformations we used was: 

O = M + Z + KxQ + P color +A,xT + P X + P y (2) 

where P co i or is a polynomial which takes into account color 
terms of first and possibly of second order; T is the time of 
the observations; P x and P v are polynomials, with order < 3, 
of the coordinates of the stars in the CCD (the x axis is in the 
West-East direction, the y axis in the North-South direction). 

All the coefficients of the transformation (0 were cal- 
culated fitting the observed magnitudes (in the instrumental 

2 available at: http://www.ls.eso.org/lasilla/sciops/2p2/E2p2M/WFI/ 
zeropoints 
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Fig. 2. Difference between cataloged V magnitudes of the stan- 
dard stars and those determined with the transformation Q. 

system) and the literature magnitudes (in the Johnson-Kron- 
Cousin system) of the standard stars. Fig. |2 shows the differ- 
ences between the V magnitudes of the standard stars from lit- 
erature and those computed with the transformation 0. The 
standard deviation of the distribution of the residuals is 0.031, 
smaller than the standard deviation (0.067) we obtained using 
the transformation Q. In the other bands the results are analo- 
gous: the standard deviations obtained with the transformation 
(0 were 0.036 in B and 0.029 in I, systematically smaller than 
those derived with transformation Q. 

In order to show the need to introduce a position polyno- 
mial term in transformation @, Fig. shows the residuals of 
the fit assuming no position dependence (upper panels) and 
with the position correction (lower panels). 

For the calibration of the B magnitudes, it was necessary 
to add to the transformations a color-dependent extinction term 
(B - V) x Q, with an em pirically determined coefficient equal 
to -0.016 (Stetson, 2005), in order to account for the higher 
extinction in this band. 

Tables lA~Tl and lA~2l in appendix|A|report all the terms and 
their coefficients used for the photometric calibration of all the 
CCDs. 

2.3. Astrometry 

The optical catalog obtained with the procedure described in 
the previous sub-sections comprised 32308 stars. We used the 
All-Sky Point Source Catalog from the Two-Micron All Sky 
Survey (2MASS) as a reference catalog to calculate the celes- 
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Fig. 4. Differential distributions of the total, spurious and cor- 
related detections of the stars in common between optical and 
2MASS catalogs. For tolerances > 1" almost all the identifica- 
tions are spurious. 
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Fig. 3. Residuals of the fit performed to calculate the transfor- 
mations for the V magnitudes in the CCDs #52 and #53 vs. the 
spatial coordinates. In the upper panels were used two transfor- 
mations without a coordinate polynomial, while in the lower 
panels the transformations included a coordinate polynomial. 

tial coordinates of the stars from their position on the CCDs. 
The 2MASS catal og was chosen for its good astrometric pre- 
cision (70-80 mas: ICutri etaill2003l) . First, we used the IRAF 
task CCXYMATCH to find a list of common sources between the 
reference catalog and our optical one, with a matching toler- 
ance of 0".5. In this procedure the celestial coordinates were 
projected in a plane with the tanx projection (a combination of 
tangent plane projection and polynomials). 

The plate solutions were computed fitting the celestial and 
pixel coordinates of the common stars with the task CCMAP. 
The transformations were then applied with the task SKYPIX to 
obtain celestial coordinates of the stars. 

The astrometric precision of the optical catalog was evalu- 
ated by matching the optical and 2MASS catalogs with differ- 
ent matching tolerances (ranging between 0".l and 2".0). We 
used a statistical procedure that allows us to sep arate the corre- 
lated and spurious identifications distributions ( Da miani et all 
I2006J) . To find the astrometric precision, we compared the dis- 
tributions of the spurious and correlated identifications in an- 
nuli with increasing radii (see Fig.|4}. We then studied the dis- 
tribution of only the correlated identifications, finding that our 
final astrometric precision was equal to 0".6 at 2cr (i.e. two 
truly correlated sources have a probability of ~95.4% of being 
identified). 

2.4. Comparison with previous catalogs 

We compared our optical catalog with the two previous cat- 
alogs of the region of NGC 6611 obtained by H93 and B99. 



The first, consisting of a total of 1026 stars with B, V, J, H 
and K photometry and limited to V< 17™, was based on ob- 
servations of the Kitt Peak National Observatory (KPNO); the 
latter is a compilation of previous catalogs including the one of 
H93, with suitable photometric transformations, and consisting 
of 2185 stars with V< 17'". 

To match these catalogs with ours, we first used a match- 
ing tolerance of 2", finding systematic differences in celestial 
positions between the common stars: -0".870 in right ascen- 
sion and +0".177 in declination for the catalog of H93, and 
-0".263 in right ascension and -0".674 in declination for the 
catalog of B99. We then used these values to correct the off- 
sets. Comparing the distributions of correlated and spurious 
detections, we estimated that the best matching tolerances were 
equal to 0".3 and 0".5 for the catalog of H93 and B99 respec- 
tively. 

We then calculated the differences in magnitude of the com- 
mon stars, with the results reported in Fig. [5] Comparing the 
median values of these differences (-0.065 and -0.081 in V 
band, 0.018 and 0.006 in B band) with the limit for accept- 
able systematic differences (taking into account anomal ous ex - 
tinction, flat-fielding errors, ecc.) of ~ 0.02 m JStetsonL Eooll) . 
there is not a good agreement in V band. As already seen, how- 
ever, the magnitudes of the standard stars are well reproduced 
by the photometric transformations used for the calibration and 
we will see later (Sect. 14. 2l that all the optical Color-Magnitude 
diagrams achieved in this work are self consistent. 

Fig- El shows the differences of the V magnitudes for the 
common stars between the catalogs of H93 and B99. There is 
not a good agreement in the V band between these two cata- 
logs. We ignored the differences in V magnitude between our 
catalog and the ones of H93 and B99, and considered our V 
magnitude to be the most reliable. 
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Fig. 5. Differences of the magnitudes in the V (upper panels) and B (lower panels) bands of the stars common to our optical 
catalog and the ones of H93 (left panels) and B99 (right panels). The dotted horizontal line marks the median magnitude difference 
between the two samples. 
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Fig. 6. Differences of the magnitudes in the B and V bands of 
the stars common to the catalog of H93 and B99. The dotted 
horizontal line marks the median magnitude shift between the 
two samples. 



3. Near-infrared, X-ray data and 
cross-identifications 

There are 30703 2MASS sources in the WFI FOV. Among 
them, we selected the NIR sources unaffected by strong pho- 
tometric uncertainties. We chose to discard the sources whose 
magnitudes were affected by instrumental artifacts, the ones 
with problems in the interpolation between the stellar profile 
and the PSF model or in the determination of the photomet- 
ric errors and the sources that are not resolved in a consistent 
fashion in all the bands 3 . The final number of 2MASS sources 
selected was then 25920. 

The X-ray data are taken from an archival 78 Ksec ob- 
servation with Chandra ACIS-I detector, with a FOV of 
16'.9 x 16'.9 i 



^inskv, 2000). Source detection was performed 
with PWDetecf , a w avelet-based source detection algorithm 
dDamiani et all Il997l) . With PWDetect we found 997 X-ray 
sources in the ACIS FOV, using a IkeV exposure map and 
a threshold limit of 5.1 <t s ^. v , corresponding to one spuri- 
ous detection on average. Fig. shows the X-ray sources of 
NGC 6611 overlaid on the WFI image corresponding to the 
ACIS FOV. This figure clearly shows that there is a crowded 
cluster of X-ray sources in the central region of NGC 66 1 1 . 

In order to obtain the multiband catalog of NGC 6611, 
we matched the common sources between the optical, near- 
infrared and X-ray catalogs. First, we matched near-infrared 
and optical detections, with a matching tolerance of 0".7, 



3 CC-flag equal to P , D, S and B; PH-qual equal to E, F and X 
or PH -qual equal to U only if the rdflag was equal to 6 JCutri et all 
120031) . 

4 available at: http://www.astropa.unipa.it/progetti_ricerca/PWDetect 
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Table 2. Results of the cross-identifications. 



Number of stars Optical detection NIR detection X-ray detection 



12561 yes no no 

16390 no yes no 

9156 yes yes no 

332 yes no yes 

33 no yes yes 

522 yes yes yes 



38995 Total of stars in the multiband catalog 

The X-ray sources without optical and 2MASS counterpart are 148. 




Fig. 7. Detail of the WFI image of Fig. Q] with the ACIS-I FOV 
(rotated square) and the 997 X-ray sources detected (north is 
up and east is on the left). The X-ray sources are crowded in 
the center of the cluster and in the North-East quadrant. 



slightly larger than the 2cr astrometric precision of the opti- 
cal catalog (0".6, see Sect. I2.3l l. The differential distributions 
of correlated and spurious identifications (Fig. O show that at 
the chosen threshold of 0".7 the correlated identifications are 
about twice the spurious ones. We have discarded the optical 
sources detected only in one band and without a NIR counter- 
part. The results of all the cross-identifications are reported in 
Table |2 with this choice and the tolerance of 0".7 we found 
9680 (522+9156 in Table EJi optical-NIR sources, with 355 
multiple identifications, 12893 (12561+332) optical sources 
without a near-infrared counterpart and 16423 (16390+33) in- 
frared sources without an optical counterpart. 

Then, we matched all the optical or NIR sources with the 
X-ray sources. We could not use a constant matching tolerance 
for all X-ray sources because the PSF in the ACIS-I images 
is not uniform: the FWHM is ~ 0".5 on-axis and becomes 
much larger near the edges of the field. We used a tolerance 
proportional to the position uncertainties derived by PWDetect 



(equal to n x cr x , where cr x is the uncertainty). To obtain the 
distribution of the spurious detections we performed 40 cross- 
identifications with different values for n and with a set of dif- 
ferent rigid translations of the X-ray sources coordinates. We 
set the value of n at 1.5 comparing the distribution of spuri- 
ous and total identifications. With this matching tolerance we 
found 887 X-ray sources with optical or near-infrared coun- 
terparts, including 78 multiple identifications. Table[3]shows a 
portion of the multiband catalog: the first column reports the ID 
of the stars; the second and third their celestial coordinates; the 
columns from 4 lh to 15 report the optical and NIR magnitudes 
with their uncertainties; the 16'' 1 is a flag indicating if the star 
has an X-ray counterpart; the 17'' 1 and 18'' 1 columns are the^ 2 , 
calculated in the optical PSF fitting procedure, and SHARP val- 
ues, a parameter calculated in the source detection procedure 
and that takes into account the sharpness of the source, allow- 
ing one to d iscard the exte nded sources and the events due to 
cosmic rays ( Stetson Ul987l) : the 19'' 1 column is a flag indicating 
if the source has a multiple identification between the catalogs 
matched to obtain the multiband catalog. 

4. The color-color and color-magnitude diagrams. 

Using the combined X-ray and optical data of the catalog pro- 
duced in this work, extended to much fainter stars (V< 23 m ) 
than previous works, we estimated new values for distance, age 
spread and average extinction of NGC 661 1. 

4.1. Distance of NGC 661 1 

Fig.|8]shows the V vs. V — I diagram of those stars in the WFI 
FOV having errors in V and V — I smaller than 0.1 m . To esti- 
mate the distance of NGC 661 1 we took advantage of the fact 
that the cluster is within the Eagle Nebula. Because of the pres- 
ence of dust associated with the nebula, the cluster stars are 
more extinguished than the foreground stars, not obscured by 
the cloud. The background main sequence stars are even more 
reddened, so in the diagram in Fig.[8]they are not visible or are 
shifted in the direction of the reddening vector. For this reason 
the main sequence stars up to the distance of the cloud in the 
diagram in Fig.[8]form a locus limited by the ZAMS (her e an d 
in the following we use the Siess isochrones: ISiess et al.l2 000>. 
We then estimated the distance of the cluster and the average 
extinction of the foregr ound stars fitting the ZA MS to this lo- 
cus, as in the study of IPrisinzano et all d2005l) of the cluster 



Table 3. Part of the electronic multiband catalog. 



ID 


RA(J2000) 


DEC(J2000) 


B 


0"B 


V 




/ 


a- 1 


J 


crj 


H 




K 




X-ray emission" 


x 1 


SHARP 


mul* 


1760 


274.847168 


-13.750009 


23.40 


0.14 


21.26 


0.04 


18.91 


0.04 


NA 


NA 


NA 


NA 


NA 


NA 


N 


1.392 


-0.028 





1761 


274.730042 


-13.749800 


22.46 


0.03 


20.10 


0.04 


17.33 


0.01 


15.42 


0.05 


14.27 


0.06 


13.74 


0.08 


N 


1.563 


0.000 





1762 


274.753571 


-13.749639 


22.60 


0.04 


20.40 


0.02 


17.53 


0.02 


15.54 


0.06 


14.58 


0.06 


13.53 


NA 


Y 


1.811 


0.601 


1 


1763 


274.733673 


-13.749538 


21.98 


0.01 


19.79 


0.02 


16.93 


0.01 


NA 


NA 


NA 


NA 


NA 


NA 


Y 


1.563 


0.119 





1764 


274.727478 


-13.749289 


23.21 


0.10 


20.91 


0.04 


17.90 


0.01 


16.25 


0.11 


14.04 


NA 


12.60 


NA 


N 


1.546 


0.059 





1765 


274.752350 


-13.749180 


NA 


NA 


23.03 


0.19 


19.75 


0.07 


NA 


NA 


NA 


NA 


NA 


NA 


N 


1.902 


-1.558 





1766 


274.828247 


-13.749168 


22.88 


0.05 


21.09 


0.03 


18.74 


0.02 


NA 


NA 


NA 


NA 


NA 


NA 


N 


1.381 


0.203 





1767 


274.731903 


-13.748953 


21.78 


0.05 


19.48 


0.01 


16.80 


0.01 


14.83 


0.11 


13.68 


0.09 


13.05 


0.06 


Y 


1.573 


-0.091 





1768 


274.714508 


-13.748638 


NA 


NA 


NA 


NA 


NA 


NA 


16.12 


0.12 


13.77 


0.06 


12.67 


0.04 


N 


0.000 


0.000 





1769 


274.827972 


-13.748659 


NA 


NA 


NA 


NA 


NA 


NA 


16.01 


0.10 


14.65 


0.06 


13.97 


0.06 


N 


0.000 


0.000 





1770 


274.737183 


-13.748060 


15.56 


0.01 


14.66 


0.01 


13.60 


0.01 


12.70 


0.02 


12.35 


0.02 


12.12 


0.02 


Y 


1.695 


0.049 





1771 


274.738922 


-13.747979 


23.97 


0.26 


22.11 


0.10 


18.67 


0.02 


16.44 


0.13 


15.13 


NA 


14.30 


NA 


N 


1.609 


0.255 





1772 


274.804474 


-13.747979 


24.32 


0.23 


22.31 


0.07 


19.55 


0.04 


NA 


NA 


NA 


NA 


NA 


NA 


N 


1.534 


0.491 





1773 


274.817169 


-13.747875 


20.79 


0.01 


19.32 


0.01 


17.60 


0.01 


16.15 


0.11 


15.27 


0.11 


14.25 


NA 


N 


1.162 


0.033 





1774 


274.829285 


-13.747687 


24.58 


0.20 


22.56 


0.30 


19.63 


0.11 


NA 


NA 


NA 


NA 


NA 


NA 


N 


1.647 


0.013 





1775 


274.726166 


-13.747549 


20.17 


0.01 


18.12 


0.01 


15.63 


0.01 


NA 


NA 


NA 


NA 


NA 


NA 


N 


1.651 


0.073 





1776 


274.722961 


-13.747337 


23.17 


0.11 


20.97 


0.03 


18.75 


0.04 


NA 


NA 


NA 


NA 


NA 


NA 


N 


1.499 


0.028 





1777 


274.759338 


-13.747359 


23.16 


0.10 


21.39 


0.04 


18.95 


0.07 


NA 


NA 


NA 


NA 


NA 


NA 


N 


1.485 


0.183 





1778 


274.756104 


-13.747284 


22.18 


0.07 


20.35 


0.01 


18.26 


0.03 


NA 


NA 


NA 


NA 


NA 


NA 


N 


1.311 


0.057 





1779 


274.720978 


-13.747194 


NA 


NA 


22.64 


0.16 


18.94 


0.03 


NA 


NA 


NA 


NA 


NA 


NA 


N 


1.816 


-0.549 





1780 


274.710754 


-13.747069 


23.03 


0.27 


21.00 


0.05 


17.75 


0.03 


NA 


NA 


NA 


NA 


NA 


NA 


Y 


1.601 


0.010 





1781 


274.806885 


-13.747075 


22.61 


0.10 


20.46 


0.02 


17.45 


0.01 


15.37 


0.06 


14.37 


0.05 


14.02 


0.08 


N 


1.630 


0.145 





1782 


274.764679 


-13.746936 


24.01 


0.29 


22.23 


0.14 


19.10 


0.08 


NA 


NA 


NA 


NA 


NA 


NA 


N 


1.530 


0.050 





1783 


274.707275 


-13.746615 


NA 


NA 


21.92 


0.37 


18.50 


0.04 


14.96 


0.05 


13.40 


0.05 


12.43 


0.04 


N 


1.712 


-0.387 





1784 


274.804321 


-13.746662 


NA 


NA 


22.59 


0.25 


20.08 


0.07 


NA 


NA 


NA 


NA 


NA 


NA 


N 


1.544 


-0.059 





1785 


274.746338 


-13.746284 


23.35 


0.20 


20.98 


0.03 


17.88 


0.01 


15.61 


0.07 


14.55 


0.05 


13.99 


0.06 


N 


1.510 


0.156 





1786 


274.791351 


-13.488843 


NA 


NA 


21.92 


0.11 


18.65 


0.05 


NA 


NA 


NA 


NA 


NA 


NA 


N 


2.398 


-0.308 





1787 


274.801392 


-13.745899 


23.08 


0.08 


20.70 


0.02 


17.75 


0.01 


15.46 


0.05 


14.31 


0.03 


13.48 


0.03 


N 


1.903 


0.086 





1788 


274.726746 


-13.745795 


24.59 


0.38 


22.33 


0.09 


19.48 


0.04 


NA 


NA 


NA 


NA 


NA 


NA 


N 


1.612 


0.256 






\fx 2 and SHARP are equal to then the source was observed in only an optical band, and was not possible to calibrate the magnitude. 
NA is for "Not Available". 



" If equal to Y then the star has an X-ray counterpart. 

* If equal to 1 the star is characterized by a multiple identification. 
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Fig. 8. V vs V - 1 diagram of the stars in the WFI FOV. Crosses 
mark the X-ray sources and identify the c luster locus. The thick 
solid line is the ZAMS from lSiess et all (EoOO) at the distance 
estimated for the cluster and with the average extinction that 
we have derived for the foreground stars (Ay = 1.45'")- The 
dashed lines are the isochrones at 0.1, 0.25, 1, 2.5, 3 and 5 
million years, with the average extinction derived for the clus- 
ter members (Ay = 2.6 m ). The extinction vecto rs are obtained 
from the relations of Munari & Carrarol {l996) and are drawn 
both for the canonical reddening law and for the one of the 
cluster (see Sect. 14. 2> . 
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Fig. 9. Histograms of Ry calculated with the U magnitudes and 
the procedure (described in the text) of H93 and our B, V and 
K magnitudes. The left panel shows the distribution of the stars 
not detected in the X-ray observation, while the right panel that 
of the stars detected. In all the panels we used all the stars com- 
mon to our catalog and to H93 and with negative Q. 

4.2. Reddening law and age spread. 

It is well known that young stars are characterized by an intense 
X-ray emission. We can assume that the large majority of the 
997 X-ray sources detected in the ACIS observ ation are young 
stars, m embers of NGC 661 1, as in the study of Da miani et alJ 
J2004 of the similar young cluster NGC 6530. 

All previous works suggest that the reddening law toward 
NGC 661 1 is slightly different to the canonical law (Ry = 3.1). 
H93 found Ry = 3.75, and this value was used in most later 
works. To determine Ry, H93 calculated the Johnson index Q 
for all the stars in their catalog w ith emission in K. Then they 
used the relations (Johnson, 1966): 



(B - V) = 0.332 x Q 



(V-K) = l.05xQ 



(3) 



which hold for main-sequence stars with negative Q (spec- 
tral class AO or ear lier); then they obtained the value of Ry 
from the equation of Steenma n & Thd ll991): 



Ry = 1.1 X Ey_ K /E B _y 



(4) 



N GC 6530. We used the ext inctions laws for the optical bands 
of lMunari & Carrarol dl996h . with Ry=3.1. 

We found that the distance of the cluster is about 1750 pc 
and the average extinction of the field stars is Ay — 1 .45"'. The 
distance of NGC 6611 we estimated is then slightly smaller 
than the values of H93 (obtained from the distance modu- 
lus of 36 OB stars) and of B99 (obtained with 72 probable 
cluster members, using a membership criterion based on the 
pr oper motions), while it is in good agreement with the value 
of iBonatto et alJ ll2006h . obtained with 2MASS data, and with 
the value reported in the Catalog of Galactic Open Clusters of 
iKharchenko et alJ J2005I) . The difference between the distance 
values found in this work and in H93 is compatible with the dif- 
ference in the V magnitude between the two catalogs. Indeed 
we obtained a distance of about 2 kpc with the procedure de- 
scribed above after correcting the V magnitudes for the median 
value of the differences shown in Fig. [3] 



and found a bimodal distribution of Ry: one peak is at the 
canonical value (Ry = 3.1, field stars) and the other at Ry = 
3.75 (cluster members). It has been shown (see Sect. 12. 4> that 
the V magnitudes obtained in our work are not consistent with 
those of H93. The K magnitudes are not in good agreement 
either, with systematic differences as for the V magnitudes. We 
chose then to repeat the procedure of H93, using our B, V and 
K magnitudes and the U magnitudes of H93. Fig.|5]shows two 
histograms of Ry, obtained for the stars without and with X- 
ray emission. The left panel clearly shows that the peak of the 
distribution of no X-ray sources is at Ry ~ 3, while for the 
X-ray emitting stars, largely dominated by cluster members, 
the peak is at Ry ~ 3.5, with a median value of Ry = 3.27, 
that we adopted in the following (right panel). This confirms 
the anomalous reddening law toward NGC 6611, even if less 
extreme than in H93. 

Once having chosen the reddening law for the cluster mem- 
bers, it was possible to redden the theoretical isochrones, from 
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ISiess et alJ j2000), in the optical Color-Magnitude diagrams 
and to estimate the age of the X-ray sources. In most of the 
previous works the age of PMS population of NGC 661 1 was 
estimated to be smaller than 1 Myr. For example, H93 found 
that the PMS stars with masses M > 3M have an age ranging 
from 0.25 to, at least, 1 Myr. As shown in Fig. [8] our data sug- 
gest that pre main-sequence stars of NGC 6611 with masses 
smaller than 1 M (with a V magnitude greater than V=17.5 m ) 
may have an age in 0. 1-3 Myr range, with a median value of 1 
Myr. 

A rough value for the average extinction for the stars of 
NGC 661 1 (A v = 2.6 m ) was estimated by fitting the bright part 
of the isochro nes older than 1 mil lion years and of the ZAMS 
from iKenvon & Hartmannl d 1995b to the bright stars with X- 
ray emission. The cluster parameters (distance, extinction and 
range of age) derived here from a V vs. V-I diagram are consis- 
tent with the parameter obtainable from the other optical Color- 
Magnitude diagrams. 

5. Stars with circumstellar disks 

We have selected the stars with circumstellar disks by searching 
for infrared excesses in their emission. To do this, we used the 
extinction-independent indices analogous to the Johnson index, 
as defined in D amiani et alJ d2Q06l) . The general definition of 
these indices is: 

Qabcd = (A-B)-(C-D)x E A - B IEc-d (5) 

where, A, B, C and D are magnitudes in four bands and 
Ea-b and Ec-d are the corresponding color excesses. The ra- 
tios between the reddening coeffi cient s was obtained from the 
extinction laws of iMuna ri & Car rarol lll996l) . for the optical 
bands, and from lRieke&Lebofskvlll985l) for the NIR colors. 
We chose these laws by comparing the extinction vectors with 
the slope of the locus of the background stars in the infrared 
Color-Color diagrams. 

We made the assumption that the V-I color is a character- 
istic of the photospheric emission of the stars. This assumption 
may be n ot strictly valid for stars char acterize d by strong ac- 
cretion (Haitig^netall[l993Enyo^& Hartmann, 1990). We 
then defined five indices to reveal excess emission in the NIR 
bands: Qvin, Qvuh, Qvijk, Qvihk, Qjhhk- 

Fig. 1101 shows the diagrams used for the selections of the 
stars with infrared excess. To be conservative, we performed 
the selection only for stars with errors, in the considered col- 
ors, smaller than 0.15'". Since the Q indices are independent of 
the extinction, the extinction vectors are parallel to the x-axis 
in Fig.[TO| We compared the indices of the stars with a lower 
limits for the photospheric emission, chose n from the normal- 
stars colors by Kenvon & Hartmann ( 1995). Only stars with an 
index smaller than the limit by more than 3 times the <tq are 
considered to show NIR excess. In order to account for spu- 
rious cross-matches between infrared and optical catalogs that 
would affect the selection of the stars with infrared excess, we 
calculated the Q indices only for the optical-NIR stars identi- 
fied within a matching radius of 0".3. 



Table 4. Number of stars in WFI and ACIS FOV and the num- 
ber of X-ray sources with excesses in the infrared bands speci- 
fied in the first column. 



Bands 


WFI FOV 


ACIS FOV 


X-ray 


J 


12 


5 





H 


7 


1 





K 


234 


138 


71 


J + H 











J + K 


5 


4 


1 


H + K 


67 


40 


19 


J + H + K 


17 


9 


3 


Total 


342 


197 


94 



Table |4] summarizes the number of stars and X-ray sources 
(Nx-rav) with excess in the NIR bands indicated in the first col- 
umn. There are very few stars with excess only in the J or H 
bands (first and second rows), and none with excess in only 
these two bands together (fourth row). The excess in K band is 
more typical of the presence of an inner disk. The nature of the 
sources with excess only in J and H might be explained with ad- 
ditional data at longer wavelength. Also only ~50% of the stars 
with a circumstellar disk are X-ray sources. The selection of 
young stars with the X-ray emission may be incomplete when 
applied to a star formation region, with a large number of stars 
with a circumstellar disk. The lack of X-ray emission in a sub- 
stantial fraction of stars with disk is consistent with the higher 
Lx of WTTS than CT TS observed in other star forming re- 
gions, su ch as in Orion (tPreibisch et al.lEo05HFlaccomio et all 
2003a b). 

6. Membership and cluster structure 

Having detected the stars with a circumstellar disk, we are able 
to identify the cluster members. The ACIS observation detected 
997 X-ray sources (see Sect.[3J, mostly members of the cluster. 
148 of the X-ray sources have neither an optical nor 2MASS 
counterpart (see Table We do not consider these stars as 
cluster members, even if their large number and their clustering 
may suggest that at least some of them are members embedded 
in the cloud. We have also excluded 18 X-ray sources without 
a circumstellar disk, because in the V vs. V-I optical dia- 
gram they have a position compatible with young foreground 
MS stars. Adding to this sample the stars with a circumstellar 
disk that are not X-ray sources, we found 1 122 cluster members 
in the WFI FOV (990 in the ACIS FOV). 

The number of cluster members found in this work is larger 
than in the previous works. H93 found a total of 174 clus- 
ter members, using a criterion that allowed them to consider 
only stars with spectral type earlier than AO. In Fig. [8] the V 
magnitude of the stars with spectral class earlier than AO is 
brighter than 15'", and X-ray sources are almost all fainter than 
this limit. Using a membership probability distribution based 
on the stars' projected position and proper motion, B99 found 
87 highly probable and 525 probable members of NGC 6611. 
In addition to the fact that our catalog is deeper than the 
one of B99 (see Sect. 12.41 the accuracy of a proper motion 



10 



M. G. Guarcello et al.: Spatial distribution of the circumstellar disks in NGC 6611 




1 2 3 4 5 0.0 0.5 1.0 1.5 2.0 2.5 

I-J J-H 




0.0 0.5 1.0 1.5 2.0 2.5 

H-K 



Fig. 10. Diagrams of the extinction-independent indices vs. the various photometric colors, used for the selection of the stars with 
NIR excess in the WFI FOV. The circles identify those stars with indices significantly smaller than the limit for the photospheric 
emission (dashed lines); the X-ray sources are represented by +, while the stars with multiple identifications with small squares. 
The thick solid lines are the photospheric stellar colors. In the Qjhhk (lowest left) diagram the thicker circles mark the stars with 
excess but without ontical emission, that are candidates embedded nrotostars. 
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a = 140 ± 8 N pc~' 
r =1-40 ±0.08 pc 



R (parsec) 

Fig. 11. Radial density profile found for the cluster members. 
The points are the observed density in circles with the radii 
shown in abscissa; the dotted line is the best-fit 2-parameters 
King's profile (equation[6j. The best-fit parameters are also in- 
dicated. The error bars are poissonian uncertainties. 

membership criterion is questionable at the large distance of 
NGC 6611. 

With the sample of cluster members we defined, we studied 
the radial density profile of the cluster. We used the 2-parameter 
density profile of King ( 1966): 



<T(f) = 



0-Q 



1 + (r/r core ) 



(6) 



where ctq is the central cluster density, r is the distance of 
the stars from the cluster center and r core the core radius. We 
defined a rough cluster center, counting the cluster members 
in squares with sides of 0.3 parsec and taking the most pop- 
ulated square as the center. The coordinates of this center are 
a = 18 : 18 : 41 and 6 = -13 : 47 : 23. Then we calculated 
the projected cluster member density (number of members per 
parsec -2 ) in concentric circles with increasing radius. To find 
the central density and core radius, we performed a non-linear 
fit between our distribution and the King's radial density pro- 
file. We found a central density (cr ) of 140 + 8 stars/parsec 2 
and a projected core radius (r core ) equal to 1.40 + 0.08 parsec 
(see Fig. II It . We repeated this procedure for the stars with a cir- 
cumstellar disk and the X-ray sources separately. We again ob- 
tained a good agreement between the observed and the King's 
profile, with a core radius of 2.34 + 0.44 parsec for the stars 
with a circumstellar disk and 1.30 + 0.08 parsec for the X-ray 
sources. This result suggests that the distribution of stars with 
a circumstellar disk is less concentrated near the center of the 
cluster, where most of OB stars are, than the distribution of the 
X-ray sources. 

The radial density pro file and the dynami cal state of the 
cluster has been studied by B onatto etal J (120061) . but using only 



o -13.80 - 




274.80 274.70 274.60 274.50 

a (J2000) 

Fig. 12. Color map of cluster members in the ACIS FOV (de- 
limited by the dotted line) in squares with sides of 0.3 pc. 

2MASS data and a membership criterion based on the NIR col- 
ors. The good agreement between the observed and the King's 
density profile suggested to Bo natto et alJ ( 120061) that the in- 
ternal region of NGC 6611 has reached some level of energy 
equipartition. Using the number of stars that they found in the 
core (120), the core radius (0.7 pc) and the definition of relax- 
ation time: 



h-elax — 771 77 ^ 

SlnN cr v 



(7) 



where N in total number of cluster members and cr,, 
the velocity dispersion, with the typical va l ue of 3 Km s _1 , 
dBinnev & Merrifieldl Il998l) . .Bonatto et all ( .20061) found the 
relaxation time to be 0.7 Myr. This value of t re [ ax is about half 
the age of the cluster that they estimated (1.3 Myr), and this 
reinforced their hypothesis about the energy equipartition level 
reached in the cluster's inner region. 

Using the value of r core found in this work, and the value 
of N for the core, estimated to be equal to 448 by our number 
of cluster memb ers and the I MF found in the study about the 
Orion cluster of iMuench et alJ (.20021) . we found a t re i ax of ~ 
4.2 Myr, greater than the upper limit of the age of the X-ray 
sources (3 Myr, see Sect. l4.2> . The hypoth esis that the core of 
the cluster has reached some level of energy equipartition is 
therefore not supported by our data. This is also confirmed by 
the significant asymmetry of the spatial distribution of cluster 
members, with very few cluster members toward the South- 
West (Fig. El). 

7. Influence of the OB stars' radiation on the 
circumstellar disks. 

We investigated the effects of the UV radiation from the nu- 
merous OB stars of NGC 66 1 1 on the evolution of the circum- 
stellar disks of the lower-mass cluster members. To do this, we 
calculated the bolometric flux emitted by all the massive clus- 
ter members earlier than B5 and incident on the positions in the 
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ACIS FOV of stars with excess in the K band. In order to take 
into account the inhomogeneous spatial distribution of the clus- 
ter, we compared the distribution of stars with disks with that 
of X-ray sources. We have not considered the absorption of the 
UV radiation of the hot stars due to the intracluster medium. 
Since the difference in extinction between the foreground stars 
and the X-ray sources is about Ay = 1 .3 and the UV radiation 
penetrates to a depth Ay ~ 1, this effect is probably impor- 
tant only for the cluster members far away from the center of 
NGC 661 1, where most of the OB stars are. We also used the 
projected distances between the stars. We will discuss later the 
effects of this approximation on our results. 

We have binned the values of incident fluxes in four ranges, 
and determined the ratio between the numbers of K-excess stars 
and X-ray sources without a circumstellar disk for each range. 
This normalization is needed to account for the spatial struc- 
ture of the cluster and it is based on the assumption that the 
structure is well traced by the X-ray sources (the richest group 
of members). The results are shown in Fig.^] where it is evi- 
dent that the fraction of K-excess stars is significantly larger in 
the first bin (i.e. at low values of incident flux) than in the other 
bins. 

If we could take into account the absorption of UV radi- 
ation by the intracluster medium and if we could use the real 
distances between the stars, instead of the projected ones, we 
would obtain lower values of incident fluxes, so the first bin 
would become more populated, reinforcing our results. Our 
analysis is therefore conservative in this respect. The spatial 
distribution of the K-excess stars is then anti correlated with 
the spatial distribution of OB stars. 

The spatial distribution of stars in the WFI FOV with ex- 
cess in the K band is shown in the left panel of Fig. [H] It is 
evident that the cluster is more extended than the ACIS FOV, 
with star-formation events in the outer regions of the nebula, 
and more dispersed than the OB stars (right panel); the spatial 
distribution of the X-ray sources is shown in Fig. [71 

7.1. Interpretation of the results 

We have found that the spatial distribution of the stars with 
a circumstellar disk in NGC 6611 is anti-correlated with the 
distribution of OB stars. We have shown that this is not due 
to energy equipartition. We have also verified that the possi- 
ble alternative interpretation, that stars close to massive stars 
are older than the remaining population (therefore with more 
evolved disks) is not supported by their position in the opti- 
cal color-magnitude diagrams, so there are no effects due to 
sequential star formation. 

We then explain the spatial distribution of stars with a cir- 
cumstellar disk with the hypothesis that UV radiation from OB 
stars causes the photoevaporation of the disks close to massive 
stars, altering the evolution of the circumstellar disks. 

Th e difference between our result and that of O liveira et alJ 
who did not find any influence of massive stars on disc 
frequency for NGC 661 1, is mainly due to the different mem- 
bers selected. Our X-ray based criterion allows us to also select 
WTTS, that are very difficult to identify with other methods. 



The process of photoevaporation works from outside 
(where the gravitational attraction of the central star is smaller) 
inward, up to the gravitational radius (r P ) , where the escap e 
velocity is equal to the sound speed (see Hollenbach 1994). 
For example, in a star with a mass of 0.2 M & with a circum- 
stellar disk heated by the incident UV radiation to 1000 K 
(FUV regime), the gravitational radius is at 20 AU from the 
central stars; in the EUV regime, when the circumstellar disk 
can be heated up to 10 4 K by the incident UV radiation, the 
gravitational radius is at 2 AU. For distances from the central 
star smaller than the gravitational radius, the photoevaporated 
gas produce s a sta ble atmosphere around the star. However, 
IWoods et alJ II 19961) have shown that if the evaporating gas 
reaches supersonic velocities in the region inside the gravita- 
tional radius, then the photoevaporation occurs even closer to 
the central stars (up to about 0.2xr g ). If the outer regions of the 
circumstellar disk are photoevaporated, accretion cannot drain 
material in the inner region replenishing the material dissipated 
by the stellar radiation or accreted on the star's surface. This 
process then may also cause the destruction of the inner region, 
responsible for the emission in NIR bands. 

We have used only the current position of the stars to cal- 
culate the flux from the massive stars incident on the circum- 
stellar disk, and we have not considered the orbit of the stars 
during the evolution of the cluster. Because of the young age 
of NGC 6611, this can be relevant only for the stars in the 
core, while the stars in the outer region have not had enough 
time to pass near the center where more OB stars are clustered. 
However, the exact computation of the incident UV radiation, 
integrated in the orbit during the lifetime of the stars with a 
circumstellar disk is beyond the scope of this paper. Using N- 
body simulations of the evol ution of clusters wit h a population 
between 100-1000 members. lAdams et alJ d2006l) claimed that 
the destruction of the circumstellar disks and young planetary 
systems are rare events. This result is not applicable in the case 
of NGC 661 1, with a larger population (1 122 members consid- 
ering only masses > O.4M ). 

8. Summary 

We have studied the young open cluster NGC 66 1 1 to derive 
the spatial distribution of the stars with circumstellar disks and 
its correlation with the population of OB cluster members. We 
found that the stars with excess in K, due to the presence of 
the circumstellar disks, are relatively more frequent at larger 
distances from OB stars. We interpret this fact as due to a fast 
disk evaporation induced by the incident UV radiation emitted 
by the numerous massive stars of NGC 661 1. 

We built a deep optical catalog compiled for the region of 
NGC 661 1, using a WFI observation in bands B, V and I, and 
comprising 32308 optical sources with magnitudes V< 23 m . 
Using our optical catalog, a CHANDRA/ACIS observation 
with 997 X-ray sources, and the 2MASS data in the WFI FOV, 
we built a multiband catalog comprising 38995 optical, NIR 
and X-ray sources. 

The catalog allowed us to improve our knowledge about 
NGC 6611. We estimated the distance of the cluster equal to 
about 1750 parsec, a value smaller than that of the work of H93 
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Fig. 13. Histogram of the ratio between the number of the stars with excess in K band in the ACIS FOV (upper numbers) and the 
X-ray sources without circumstellar disk (lower numbers) vs. the incident flux from OB stars. We also show the error bars of the 
ratio. 



K-excess stars OB stars 




Fig. 14. Spatial distribution of the stars with excess in the K band in the WFI FOV (left panel) and of the OB stars members of 
NGC 661 1 (right panel). The rotated square marks the ACIS FOV. 



and B99, but similar to the one found in more recent works; we 
also confirmed the presence of an anomalous reddening law in 
the direction of NGC 661 1, with R v = 3.27. Using the X-ray 
emission as a membership criterion, we found a large number 
of cluster members with masses < 1 M Q . We found stars with 
ages up to 3 Myr together with very young stars with ages of 
about 0. 1 Myr. Most of the previous works identified only stars 
with ages less than 1 Myr. 

We identified 342 sources with excesses in the NIR bands 
in the WFI FOV using suitable color indices, independent of 
the extinction. In the ACIS FOV, 47% of the stars with circum- 



stellar disks have an X-ray counterpart (94 X-ray sources with 
192 stars with K band excess in the ACIS FOV). The stars with 
circumstellar disks and those that are X-ray sources form the 
most numerous and the deepest list of members of NGC 661 1 
(1 122 cluster members in total). 

Using this member list we have studied the radial density 
profile of the cluster, finding a core radius of 1 .40 + 0.08 parsec. 
The spatial distribution of the stars with a circumstellar disk 
is less concentrated near the cluster center, where most of the 
OB stars are, than the X-ray sources. The fraction of stars with 
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disks increases with decreasing incident UV flux from early 
type stars. 

There are 194 stars with a circumstellar disk outside the 
ACIS FOV, showing that the Eagle Nebula is characterized by 
star-formation events. 
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Appendix A: Coefficients of the solution for the photometric calibration 



Table A.l. Zero-points, coefficients of the colors and temporal terms for the CCDs, used in the transformation (fJJ- 
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A„ is the coefficient of the color term with exponent equal to n. 
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Table A.2. Coefficients of the position's polynomials used for the photometric calibration in all the CCDs. 
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